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O 1 The solution to the non-forward BFKL equation in the Leading Logarithmic approxima- 

tion is expressed in terms of a sum of iterations of its kernel directly in transverse momentum 
and rapidity space. Several studies of the non-forward solution are performed both at the 
level of the gluon Green's function and for a toy cross-section, including an analysis of the 
diffusion properties as found in this approach. The method developed in this paper allows 
for a direct inspection of the momenta in the BFKL ladder, and can be applied to solving 
the non-forward BFKL equation to next-to-leading logarithmic accuracy, when the corre- 
sponding kernel is available. 
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1 Introduction 



An interesting framework to study the behaviour of QCD scattering amplitudes in the limit 
of large centre-of-mass energies y/s and fixed momentum transfer yf—i is the Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) formalism [1]. The BFKL framework can be applied to processes 
characterised by a colour octet exchange as well as colour singlet exchange (diffractive processes) . 
The case of t = colour singlet scattering corresponds to forward scattering, whereas the case 
of t ^ is called non-forward scattering. The optical theorem relates the amplitude for forward 
diffractive scattering to the amplitude for a colour octet exchange, which will be exploited in the 
check of some of the results derived in this paper. If the momentum transfer is perturbative, i.e. 
—t 3> Aq CD , it is possible to use the non-forward BFKL equation to study high-t diffraction in 
the high energy limit, which is characterised by a final state with two systems with transverse 
momentum —t and far apart in rapidity. The colour singlet exchange in the non-forward case 
results in a rapidity gap in jet activity. The non-forward BFKL equation thus provides a useful 
theoretical framework to study diffractive physics from first principles in QCD. 
The differential cross-section for diffractive scattering, 

da \A(s,t)f 

tit 16tts 2 ' 1 j 

can be described within the BFKL framework by a factorised scattering amplitude of the form 

\Ms,t)\ 

s 

where Y is the rapidity separation of the scattered probes, <3?A(k a ,q) and 3>B(k&,q) are the 
process-dependent impact factors and the four-point gluon Green's function, /(k a , k^, q, Y), is 
universal. For illustration, in Fig. ^ a typical dominant contribution to this colour singlet ex- 
change in the high energy limit is shown. As explained below in Sec. [2] the non-forward BFKL 
equation describes the evolution of /(k a ,k&,q, Y) as a function of the rapidity separation Y. 
This formalism can be applied, for example, to the study of the diffractive production of vector 
mesons in photon-proton collisions at HERA for large center-of-mass energies and transverse 
momentum squared \t\ > Aq CD . In this kinematical region the Leading Logarithm (LL) terms 
(a s ]ns/\t\) n generated in the perturbative series must be resummed. Non-perturbative contri- 
butions are included in the proton parton densities and in the meson light-cone wave function 
present in the corresponding impact factor. For the p, cj) and J/ip mesons the transverse mo- 
mentum spectrum and spin density matrix elements have recently been measured at HERA [2]. 
From the theoretical side, there has been an intense activity in the study of these processes, 
see e.g. Ref. [3]. The case of photon dissociating to a photon has also been studied in Ref. [4]. 
Another example of the application of this approach is the description of events with interjet 
rapidity gaps in photon-hadron and hadron-hadron collisions [5]. 

The original analytic solution to the non-forward LL BFKL equation [6] proved to be sig- 
nificantly more complicated than its forward counterpart. In order to exhibit the conformal 
invariance that proved vital in solving the forward BFKL equation, it is necessary to perform a 
Fourier transform from transverse momentum space to impact parameter representation. The 
non-forward BFKL equation at next-to-leading logarithmic (NLL) accuracy (a s (a s lns/\t\) n 
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Figure 1: A Feynman diagram contributing to the LL approximation for quark-quark scattering 
with colour singlet exchange. 

terms) will be significantly more complicated, with the conformal invariance being explicitly bro- 
ken by the running of the coupling (as it happens in the NLL forward case [7]). It is therefore 
important to investigate new strategies to find the solution to this equation. 

In this paper we present a new approach to obtaining the solution to the non-forward 
BFKL equation to LL accuracy. The method of solution is based on the separation of different 
contributions to the BFKL kernel by a phase space slice. Once this separation is performed it 
is possible to apply an iterative approach similar to the one presented in Ref. [8] for the case 
of the BFKL equation describing a colour octet exchange. This approach has recently been 
generalised [9] to solve this BFKL equation in the NLL approximation both in QCD and N = 4 
supersymmetric Yang-Mills theory [10] (reviews can be found in [11]). 

In the present work we concentrate on extending this iterative approach to the solution 
of the BFKL equation describing diffractive processes. The presented method of solving the 
BFKL equation directly in transverse momentum space has the benefit that it allows for a 
direct inspection of all involved momenta. In particular it is possible to study the diffusion of 
the transverse scales along the evolution in rapidity. The same method of solution can be applied 
to the non-forward BFKL equation in the NLL approximation. The corresponding analysis will 
be performed when the calculation of the non-forward NLL kernel is completed [12]. 

The structure of this paper is as follows. In section |2] the LL BFKL equation in the non- 
forward case is presented and the phase space slicing regularisation procedure is performed. In 
section |31 the non-forward LL BFKL equation is solved by iterating the kernel, obtaining in this 
way an explicit expression for the four-point gluon Green's function. In section0]the behaviour of 
the four-point gluon Green's function is analysed for different values of the momentum transfer. 
A study of the diffusion properties is performed of not just the four-point gluon Green's function, 
but also, in section |SJ for cross-sections calculated with toy impact factors. The conclusions are 
presented at the end. 
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2 The Non-Forward LL BFKL Equation 



The LL BFKL equation describing colour singlet exchange with a non-zero momentum transfer 
squared, —t = q 2 , was originally calculated in Ref. [1]. The starting point in the present work 
will be Eq. (4.16) of Ref. [13] where the integral equation for the Mellin transform in rapidity 
of the four-point gluon Green's function, f u (k a , kt, q), is given by 



wfu (ka,k 6 ,q) 



5^ (k„ - k 6 ) + g / d 2 k' 



+ 
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where k a and k^ describe the two-dimensional transverse momenta of the exchanged gluons 
in the t-channel (see Fig. ^) and we define a s = a s N c /n. The driving term 5^ (k a — k&) 
corresponds to a simple two gluon exchange. It is more convenient to rearrange the terms in 
this expression and, as the integration variable, to use the transverse momenta of the s-channel 
gluons, k = k' — k a , i.e. 



ufu (k a ,k 6 ,q) = <5 (2) (k a -k 6 



a 



(kg-q) 2 (k + kq) 2 X 

(k + k a -q) 2 k 2 
k 2 
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f u (k + k a ,k ,q) 



(k a 



k 2 I (k + k a ) 2 +k 2 (k + ka 



k 2 



f u (k a ,k 6 ,q) 



(4) 



At this stage it is possible to separate the exchange terms form the s-channel contributions by 
introducing a phase space slicing parameter A (at NLL it is useful to use dimensional regular- 
isation to show the cancellation of infrared divergences, at LL this is not needed). It is also 
convenient to use the following approximation 



/ w (k + k„,k b ,q) = / w (k + k a ,k 6 ,q) (#(k 2 -A 2 ) +9(X 2 -k 2 )) 

~ fu (k + k a , k ft , q) 6 (k 2 - A 2 ) + fu (k a , k 6 , q) 9 (A 2 - k 2 ) 



(5) 



In all the results presented in this paper we have made sure this approximation is valid by 
checking that the four-point Green's function is insensitive to the value of the slicing parameter 
for small values of A. 

Therefore, the non-forward BFKL equation for the Green's function can be written in a very 
simple form: 



(w-a;o(k a ,q,A))/ w (k a ,k 6 ,q) = <5 (2) (k a - k 6 ) 

r 

+ J ^0 (k 2 - A 2 ) £ (k a , k, q) fu (k + k a , k 6j q) . 



(6) 
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In order to write the equation in this way the following notation has been introduced: 



(k + k a -q) 2 k 2 



(7) 



and 



o, (k a ,q,A) = ^ ^ 



(A 2 -k 2 ) 1 + 



(k a -q)^(k + k a f -q 2 k 2 



(k + k a -q) 2 k 2 



(k a " q) 



(k + k a ) 2 + k 2 (k + k a -q) 2 +k 2 



(8) 



The latter expression corresponds to the Regge trajectory in our regularisation. In the case of 
q = 0, the trajectory for the colour-octet exchange should be obtained, and indeed we find 



ojq (k a , 0, A) 



a s f <rk 

W k 2 " 



;a 2 - k 2 ) 



(k+k a ) 2 + k 2 



a s ln^f. (9) 



Moreover, the whole solution to the non-forward BFKL equation has the correct limit for q — * 
as found in Ref . [8] . 

The expression for the non-forward LL Regge trajectory can be simplified if we make use of 
the forward limit in Eq. ©, i.e. 



uj (k a ,q, A) = - (oj (ka, 0, A) + lu (k a - q, 0, A)) 

« s /-d 2 k 2 |2 , ^ (k a -q) 2 (k + k a ) 2 -q 2 k 2 ' 
+ T.j k 2 "^ (k + k a -q) 2 k 2 \ 



(10) 



The last integral is negligible when A is small, and therefore the non-forward trajectory takes 
the simple form 



uj (k a ,q, A) 



(ln^+ln (ka ~ q) 
2 A 2 A 2 



(11) 



With these conventions we proceed in the next section to iterate Eq. © to find the solution for 
the non-forward four-point gluon Green's function. 



3 Solution to the Equation 

The non-forward LL BFKL Green's function is the solution to the integral equation 

1 



f w (k Q ,k fe ,q) 



lo-uo (k a ,q, A) 
d 2 k 



{s™ (ka-kfc) 



+ J ^0 (k 2 - A 2 ) £ (ka, k, q) / w (k + ka, k 6 , q) j . (12) 



5 



for A — ► 0. If this expression is iterated, the Green's function can be expressed in terms of a 
kernel per iteration acting on the initial condition with a series of poles in the uj complex plane: 



/a, (k a , k 6 , q) = 



+ 



+ 



+ • 



£ (2) (kg - kfc) 
lo - uj (k a ,q, A) 

d 2 kj 6 (kj - A 2 ) C (K, ki, q) (kg + ki - kft) 

7rkf w - w (k a ,q, A) w - a; (k a + ki, q, A) 

d 2 ki /" <i 2 k 2 6 (kf - A 2 ) e(ka,ki,q)e(k2-A 2 ) £ (k a + ki, k 2 , q) 



vrk 2 



vrk 2 



w - (j (k a ,q, A) 
5^ (kg + ki + k 2 - kg) 
' w - w (k a + ki + k 2 , q, A) 



w - ojq (k a + ki,q, A) 



(13) 



Each action of the kernel corresponds to an interaction between the reggeised gluons exchanged 
in the t-channel, building up, in this way, the LL BFKL ladder. The poles are integrated over 
when going from the w-plane to rapidity space using the Mellin transform 



/(k a ,k fe ,q,Y) 



1 

2vri 



a+ioo 



dw e" Y f u (k a ,k 6 ,q) , 



(14) 



where Y is the rapidity span of the BFKL ladder. 

This integration can be performed to finally obtain the solution of the non-forward LL BFKL 
equation for the four-point Green's function as (with yo = Y) 



/(k ffi ,k 6 ,q,Y) 
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This solution has the correct forward limit of Ref. [8] when the momentum transfer tends to 
zero. We would like to stress again that this method of solution is directly applicable to the non- 
forward BFKL equation also at NLL accuracy, with the difference that the introduction of the 
phase space slice is more conveniently performed in the language of dimensional regularisation. 



4 Analysis of the Gluon Green's Function 

To study the dependence of the non-forward LL BFKL four-point gluon Green's function on 
the transverse momentum scales we choose to integrate over all external angles and define the 
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Figure 2: The angular integrated Green's function for fixed | kit, I = 5 GeV as a function of |k a | 
for different values of the momentum transfer. The rapidity span is low, Y = 1. 

quantity 

f'I'K flit 

/(|k a |,|k 6 |,|q|,Y) = / d6 qa / d0g 6 /(k o ,k 6 ,q,Y) J (16) 

Jo Jo 

where 6 q i is the angle between the vectors kj and q. 

As a first analysis, in Fig. |2] the value of |k&| is fixed to 5 GeV and the dependence on |k Q | 
is studied. When the value of the modulus of both momenta coincides, the angular integrated 
Green's function shows a 5-functional behaviour corresponding to the two gluon exchange limit. 
This dependence is caused by the driving term of the integral equation whose influence is stronger 
for lower energies. When Y is increased from 1 in Fig. |2 to Y = 3 in Fig. |3] the influence of 
the driving term diminishes as a consequence of a larger number of effective rungs in the BFKL 
ladder. It is interesting to note that the influence of the momentum transfer q is larger in regions 
of low scales of k a . The general trend is that the four-gluon Green's function decreases with 
increasing q. Comparing Fig. |21 with Fig. it can be seen that this effect is increasing with 
rapidity. 

Let us now proceed to the study of the diffusion of the transverse scales in the BFKL ladder 
with increasing rapidity span. We choose in this analysis to concentrate on the evolution of the 
transverse scale along the left hand side of the ladder depicted on Fig.Q but, as previously men- 
tioned, one of the benefits of having the solution to the non-forward BFKL equation expressed 
in terms of explicit phase space integrals as in Eq. (|15[) is that it is possible to study any of the 
momenta in the ladder. 
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Figure 3: The angular integrated Green's function for fixed |k&| =5 GeV as a function of |k a 
for different values of the momentum transfer q. The rapidity span is Y = 3. 



Diffusion is normally studied in terms of the mean of the transverse momentum along the 
ladder, with the width of the distribution indicated by the standard deviation (see e.g. Ref. [13, 
14]). This choice is useful, since it allows analytic studies. However, a plot of the mean internal 
transverse momentum plus/minus the standard deviation as a function of the rapidity along 
the BFKL ladder fails, by construction, to display the different behaviour of diffusion to low 
and high scales. Given the possibilities offered by the solution written in terms of explicit 
phase space integrals, we therefore choose to study diffusion in terms of the average value (r) 
of r = ln((k a + h) 2 /GeV 2 ) as a function of the rapidity Y' along the ladder. Specifically, for 
a given value of k a ,kb, and Y we solve the non-forward BFKL equation according to Eq. (|15|) 
by a MonteCarlo integration method. For each configuration point in n— momenta phase space 
{ki,Ui} we can trace the evolution of r along the ladder and, at the same time, calculate the 
weight of this configuration to the total solution. In this way it is possible to calculate both the 
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Figure 4: The diffusion properties of the gluon Green's function in terms of the lines formed 
by (r), (r) + a%, and (r) — 02 along the BFKL ladder. Shown for k a = 5 GeV, k b = 4 GeV, 
<? = GeV (left) and q = 2 GeV (right), for rapidity spans of Y = 1, 2, 3. 

average value of r along the ladder, (r)(Y'), and the quantities 

POO 

2/ dT(T-(T)(Y')) 2 f(k a ,k b ,q,Y) 



(17) 



of(Y')- J< ^ Y,) 



00 



dTf(k a ,k b ,q,Y) 







•<T>(Y') 

2/ dr(r-(r)(Y')) 2 />a,^,g,Y) 



00 



dTf(k a ,k b ,q,Y) 







For a gluon Green's function /(k a ,k&,q, Y) that is symmetric in t, the lines of (r)(Y'), (t)(y') + 
o"i(Y'), and (t)(y') — cr 2 (Y') would reproduce the plot of the mean plus/minus the standard 
deviation. This is true for q = GeV, as shown in Fig. |I] Here we have plotted the above 
mentioned three lines for k a = 5 GeV, k b = 4 GeV, q = GeV, a s = 0.23, and Y = 1, 2, 3. Y' is 
rescaled to lie between and 1, so as to plot all the three cases on the same figure. 

However, for q > GeV the distribution of r is no longer symmetric, as can be seen on the 
right hand side plot in Fig. |1J where we have plotted the same quantities for q = 2 GeV. It is 
apparent that there is less diffusion to smaller scales, the average is unchanged, and the diffusion 
to scales larger than the average is almost unchanged compared to the case q = GeV. 

It is also apparent that for the set of parameters investigated here, there is no significant 
diffusion into regions where the coupling is expected to become unperturbatively large. As it is 
well known, the influence of softer scales is larger at higher energies, but in this work we confirm 
the fact that the diffusion into the infrared is drastically reduced when there is some momentum 
transfer, acting, in this way, as an efficient infrared cut-off. It will be interesting to investigate 
if this picture holds at NLL. 
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Figure 5: The impact factor &A{k a )/ka for q = GeV and m = 3.1 GeV. 



5 Analysis of a Toy Cross— Section 

Let us now turn to the study of diffractive cross-sections. To this end we need to define some 
suitable impact factors. We choose a generic example from Ref. [13] modelling the photo- 
production of a vector meson: 



<£> A (k a ,q) = a s h 2 / dpdrp(l- p) 



(q 2 p 2 -r(l — t) + m 2 ) ((k a — pq) 2 r(l — t) + m 2 



> (18) 



and similarly for <3?b- Here, h is a normalisation constant which we choose arbitrarily such that 
a s h 2 = 1. This choice obviously means that the normalisation of the toy cross-section reported 
in this paper is completely arbitrary. We have plotted this impact factor divided by k\ in 
Fig. El for q = GeV and the mass of the meson m = 3.1 GeV. One further complication arises 
compared to studies of the four-point gluon Green's function due to the integration over k a and 
kb when calculating the differential cross-section as in Eq. I[T|). The approximations of Eqs. (JHJ 
and (jllj) are valid only when A <C k a . When k a is integrated over it is therefore not possible to 
have A fixed. This situation is similar to the one encountered in Ref. [15]. This problem can be 
solved in several different ways. In the present analysis we choose to always have A < 20 k a and 
A < 1 GeV. We have checked that the results here presented do not depend on these choices. 

On Fig. |13 we have plotted the resulting cross-section of Eq. (Q) as a function of the rapidity 
span of the BFKL ladder for q = 0, 1, 2 GeV, a s = 0.2 and m = 3.1 GeV. The exponential rise of 
the cross-section as a function of rapidity is evident for all q. In the previous section we showed 
how the Green's function diminishes as the momentum transfer rises, this translates here into 
smaller cross-sections. 

In Fig. [7| we have plotted the q— dependence of the cross-section for rapidity spans of Y = 
0,1,2,5. We see an exponential fall-off with q for all Y. Furthermore, we observe that da/dt 
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Figure 6: The toy cross-section as a function of the rapidity separation of the produced vector 
mesons (m = 3.1 GeV), each produced with transverse momentum q = 0,1,2 GeV. 

at t = GeV 2 increases with Y. The curves in Fig. [7| are presented on a linear (right) and a 
logarithmic (left) scale. 

These results could have also been obtained using an analytic approach to solving the non- 
forward BFKL equation at LL accuracy. Nevertheless it is when discussing diffusion properties 
that the method presented in this paper has advantages. For the analysis of diffusion in cross- 
sections, a second issue arises compared to the study of diffusion of the Green's function. The 
perturbative scale is now set by both the meson mass m and the momentum transfer squared 
—t, while the scales k a and are no longer fixed. It therefore becomes interesting to study not 
only the distribution of the average transverse scale (r) along the ladder, but also at the ends of 
it, where it will describe the average scale of the transverse momentum connecting the impact 
factors to the gluon Green's function. 

In Fig.|H]we have plotted the distribution of the average momentum scale in a similar way to 
the one used in Fig.0]for the Green's function. It is interesting first to note that the spread along 
the chain is not significantly larger than that at the ends of the chain (for the rapidity spans 
considered here). Secondly, it is comforting to see that the typical scales remain perturbative 
for all values of the momentum transfer q, for the chosen value of the mass of the vector meson 
(m = 3.1 GeV). The average value of r for q = GeV ((r) « 3.9) corresponds to a scale of the 
internal momenta of the BFKL exchange of roughly 7 GeV. The logarithmic scale of momenta 
on Fig. ^ and Fig. |H] emphasises the region of soft momenta. It might therefore be helpful to 
study the effect of the increase in scattering momenta on the internal scales directly: an increase 
of q from GeV to 2 GeV leads to and increase in the upper (UV) lines corresponding to an 
increase in the internal momenta of roughly 2.4 GeV, while the increase in the lower (IR) lines 
corresponds to an increase of roughly 1.7 GeV. 
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Figure 7: The toy cross-section as a function of the transverse momentum of the produced 
vector mesons, for rapidity separations of Y = 0, 1, 2, 5 units of rapidity. 
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Figure 8: The range of relevant internal transverse momentum scales in terms of (r), (r) + 
01, (r) — 02 for the diffractive cross- section for a s = 0.2, m = 3.1 GeV, Y = 2,3, and q = 
0, 1, 2 GeV. Y' is the rescaled (to unity) rapidity along the BFKL ladder. 
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6 Conclusions 

We have presented a new solution to the non-forward BFKL equation at leading logarithmic 
accuracy that allows a study of diffusion properties directly in momentum space. We have 
investigated the behaviour of the gluon Green's function as a function of transverse scales, 
including a study of the IR/UV diffusion. Then we extended this study to the analysis of a toy 
cross-section. 

The presented framework is very efficient for solving BFKL evolution equations and the 
solution allows immediate insight into the momentum configurations of the evolution. We hope 
to extend the iterative method for solving the non-forward BFKL equation to next-to-leading 
logarithmic accuracy, once the appropriate integral kernel is calculated. 
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